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In this talk, we review our recent calculation for the pp — > ppir reaction near threshold in pionless 
effective field theory with a di-baryon and external pions. 
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I. INTRODUCTION 

The study of neutral pion production in proton-proton 
collision near threshold, pp — > ppir , has been inspired 
by precise measurements of the near-threshold cross sec- 
tion^. Surprisingly, the measured cross section turned 
out to be ~ 5 times larger than the early theoretical pre- 
dictions^. Subsequently, some mechanisms to account for 
the threshold experimental data have been suggested in 
model calculations^. 

Heavy-baryon chiral perturbation theory (HB^PT) is 
a low-energy effective field theory (EFT) of QCD and 
provides us a systematic perturbation scheme in terms of 
Q/A x where Q denotes small external momentum and/or 
symmetry breaking term m w and A x denotes the chiral 
scale A x = Airf^ ~ 1 GeV: f n is the pion decay constant. 
Though many works on the pp — > ppir 
threshold in HBxPT have been done^&i^ 
recent review, see Ref.— and references therein), some 
issues in theoretically describing the process have not 
been fully clarified. In the next-to-leading order (NLO) 
HBxPT calculations^^, a significant enhancement of 
the off-shell irirNN vertex function obtained from the 
NLO HB^PT Lagrangian is found. However, the two- 
body (one-pion-exchange) matrix element with the off- 
shell irirNN vertex is almost exactly canceled with the 
one-body matrix element. Thus the experimental data 
cannot be reproduced in the NLO calculations. In the 
next-to-next-to leading order (NNLO) HBvPT calcula- 
tions^^, a significant contribution comes out of the 
NNLO corrections and a moderate agreement with the 
experimental data is obtained—. However, the chiral se- 
ries based on the standard Weinberg's counting rules^ 
shows poor convergence. 

A modification of the original Weinberg's counting 
rules to account for the large momentum transfer, k ~ 
y/m^mN where and m^r arc the pion and nucleon 
masses, respectively, is discussed in Ref.»S. The produc- 
tion operators at NLO using the modified counting rules 
are estimated, and it was reported that the NLO con- 
tributions exactly cancel among themselves^. Recently, 
some detailed issues for the loop calculations, such as a 
concept of reducibility 1 -, a representation invariance of 
the chiral fields among the loop diagrams, and a proper 
choice of the heavy-nucleon propagator—, were also stud- 
ied. 

Meanwhile, it is known that the energy dependence of 



the experimental data can be well reproduced in terms 
of the final state interaction and the phase spaced. A 
"minimal" formalism to take account of these two fea- 
tures would be a pionless theory, in which virtual pions 
exchanged between the two nucleons are integrated out; 
in this pionless theory, the one-pion exchange, two-pion 
exchange and contact terms in HB%PT are subsumed in 
a contact term. Furthermore, after taking these two fea- 
tures into account, the difference between the theory and 
experiment appears in the overall factor and the exper- 
imental data can be easily reproduced by fitting an un- 
known constant that appears in a contact vertex. In this 
worki£ we employ a pionless EFT with a di-baryo n 19 i 2030 
and external pions 2 ^ to calculate the total cross section 
of the pp — ► ppir process. 



II. PIONLESS EFFECTIVE LAGRANGIAN 

An effective Lagrangian without virtual pions and with 
a di-baryon and external pions for describing the pp — > 
ppir reaction may read 
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where Cm is the standard one-nucleon Lagrangian in 
heavy-baryon formalism where the "external" pions are 
nonlinearly realized. C s is for the 1 <S'o channel di-baryon 
field, Cnstt represents the contact interaction of an exter- 
nal pion, a di-baryon and two nucleons, and C NN is for 
the two-nucleon 3 Pq channel. The effective Lagrangian 
for the two-nucleon part may rea d 19 ' 20 i 24 ' 25 
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with O 

-^<J2<JiT2T a . s a is the 1 S'o channel di-baryon field and 
o~ s is the sign factor o~ s = ±1 which we fix below, is a 



FIG. 1: Diagrams for the dressed di-baryon propagator in- 
cluding the Coulomb interaction. 
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FIG. 3: Diagrams for the P-wave NN scattering. 




FIG. 2: Diagram for the .S'-wave pp scattering amplitude with 
the Coulomb interaction. 




FIG. 4: Diagrams for pp — » ppn near threshold with the 
strong and Coulomb final state interactions and without the 
initial state interaction. 



velocity vector = (1,0) and D M is the covariant deriva- 
tive. S s is the mass difference between the di-baryon mass 
m s and two-nucleon mass, m s — 2m^ + 5 s . y s is the cou- 
pling constant of the di-baryon and two-nucleon interac- 
tion, iv = -n=T2T a <72- is an unknown low en- 

1 Vo 

ergy constant (LEC) of the (external) pion-(spin singlet) 
dibaryon-nucleon-nucleon (irsNN) interaction, ro is the 
effective range in the 1 Sq (pp) channel, and A M = ?fA%. 
C*2 is the LEC for the NN scattering in the 3 Pq channel. 

Now we calculate the S- and P-wave NN scattering 
amplitudes to fix the LECs in the two-nucleon part. In 
Fig. [TJ diagrams for the dressed 1 S'o channel di-baryon 
propagator are shown where the two-nucleon bubble dia- 
grams including the Coulomb interaction are summed up 
to the infinite order. In Fig. [21 a diagram of the S'-wave 
pp scattering amplitude with the Coulomb interaction is 
shown and thus we have the S'-wave scattering amplitude 
as 



iA s = {-iy s ^o) [W s (p)] (-iy s ip ) 
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where ipo = C^e lCT ° and Co is the S-wave Coulomb phase 
shift cto = argT(l + irj). D s (p) is the dressed di-baryon 
propagator and 5f is the renormalized mass difference. 
h(rf) = Re ip(iij) - lnr?, Reip{ij) = i] 2 YlT=i v {J+r?) ~ 7 ' 
7 = 0.5772 ■ • and 
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The S-wave amplitude A s is given in terms of the effec- 
tive range parameters as 
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where ac is the scattering length, r is the effective range, 
and the ellipsis represents the higher order corrections. 
Now it is easy to match the LECs with the effective range 



parameters. Thus we have a s = — 1 and 
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We note that the sign of the LEC y s cannot be deter- 
mined by the effective range parameters. 

In Fig. [31 diagrams for the P-wave NN scattering are 
shown. Because the momenta of the two protons are 
quite large for the pion production reaction, the two- 
proton bubble diagrams are summed up to the infinite 
order. The scattering amplitude for the 3 Po channel is 
obtained as 
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The LEC C® is fixed by the phase shift of the 3 Po channel 
at pion production threshold, 5 v (pth) — —7.5° at pth — 
^/m^mpf. Thus we have 
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III. 
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In Figs. [4] and [5j we show diagrams for pp — ► ppn near 
threshold. In diagram (a) in Fig. Q] and (c) in Fig. [51 the 
pion is emitted from the one-body ttNN vertex. In the 
diagram (b) in Fig. [Hand (d) in Fig. [51 the pion is emitted 
from the nsNN contact vertex which is proportional to 
the LEC d^> . The one-body amplitude from the (a) and 
(c) diagrams and the two-body (contact) amplitude from 
the (b) and (d) diagrams are obtained as 



+<■■) 



^9A 



4-TT L P 



3 



(c) (d) 

FIG. 5: Diagrams for pp — + pprv with the strong initial and 
the strong and Coulomb final state interactions. 




FIG. 6: Diagram for a one-pion exchange contribution to the 
pp — > pp-K process for estimating the LEC di 2 ' in the contact 
vertex. 
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FIG. 7: Estimated total cross section of pp — > ppiv as a 
function of — \q\ m .ax/m 7r . See the text for details. 
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where 2p and 2pl are the relative three momenta between 
incoming and outgoing two protons, respectively; p = |p| 
and p' = \p*\. r/' = arnjv/(2|> / ) and oj q is the energy 
of outgoing pion, uj q = \J q 2 + m 2 : q is the outgoing 
pion momentum. We note that there remain no unknown 
parameters in the amplitudes except for the LEC d\ in 
the two-body (contact) amplitude in Eq. (fl"3"|) . 

Now we estimate an order of magnitude of the LEC 
d\ from HBxPT. We here consider a one-pion-exchange 
(OPE) diagram shown in Fig. O This diagram is the 
lowest order OPE contribution in the standard Weinberg 
counting rules. We include a higher order (relativistic) 
correction to the ttttNN vertex which is found to be im- 
portant^ and is, in the modified counting rules, of the 
same order as the lowest order diagram. 

The effective chiral Lagrangian to calculate the 
isoscalar ttttNN interaction in the diagram in Fig. [51 
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where 
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and x is the relativistic correction to the term propor- 
tional to (v ■ A) 2 in Eq. (TH). The values of the LECs ci, 
C2 and C3 are fixed in the tree-level calculations 29 as 
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c 1 = -0.64, c 2 = 1.79, 



C3 



-3.90 [GeV- 1 ] . (15) 
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The value of the LEC eft from the loop diagram in 
Fig. [6J is obtained as 
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where the different signs for d\ have been obtained be- 
cause of the LEC y s in Eq. (J5J . 
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IV. NUMERICAL RESULTS AND SUMMARY 



The total cross section of pp 
calculated using the formula 



ppir near threshold is 
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p' = W\ =i ^Jm N {T- y/ml + q 2 )-q 2 Jl, and 
. T is the initial total energy T ~ 



T 2. 



p^/rriN and ~ 2p/mjv- We have expanded the 
proton energies in the phase factor in terms of I/uin 
and kept up to the 1/mjv order. .4 is the amplitude 
.4 = -4( Q+C ) + -4(6 + d) where -4( Q+C ) and „4( 6+ci ) are ob- 
tained in Eqs. (fT!?]) and |T3]), respectively. 

In Fig. [7] we plot our results for the total cross section 
as a function of r) n = q max /m^. The solid curve and long- 

dashed curve have been obtained by using d\ = ±0.140 
fm 5 / 2 fixed from the one-pion exchange diagram in Fig. [6j 
in the previous section. The LEC d^ is also fixed by 
using the experimental data as 
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where we have two values of with different signs 
The short-dashed curve is obtained by using ^ 2 ^ ltted = 



4 



—0.12 fm 5 / 2 . The dotted line corresponds to the case 
where only the contribution from the one-body amplitude 
A( a +c) is considered. The experimental data are also 
included in the figure. 

We find that the experimental data are reproduced 
reasonably well with the value of — —0.14 fm 5//2 . 
By contrast, we obtain almost vanishing total cross sec- 
tions with the value — +0.140 fm 5//2 because the 
two-body amplitude with Ji 2 ' = +0.140 fm 5 / 2 is almost 
canceled with the amplitude from the one-body contri- 
bution. On the other hand, for the whole energy range 
the experimental near-threshold cross section data are 
well reproduced with the use of the fitted parameter 
fitted _ _Q i2 fm 5/2 . We also find that approxi- 
mately a half of the observed cross section comes from 
the one-body (IA) amplitude in the pionless theory. 

In this work we calculated the total cross section for 
pp — > ppir near threshold in pionless EFT with the di- 
baryon and external pion fields. The leading one-body 
amplitude and subleading contact amplitude were ob- 
tained including the strong initial state interaction and 
the strong and Coulomb final-state interactions. After 
we fix the LECs for the NN scatterings, there remains 
only one unknown constant, eft- , in the amplitude. We 
estimated it from the one-pion exchange diagram in the 
pionful theory. Although this method does not allow us 

~(2) 

to fix the sign of d\\ we have found that one of the 

~(2) 

two choices for d^ leads to the cross sections that agree 



with the experimental data reasonably well. On the other 
hand, the whole range of the experimental data near 
threshold can be reproduced by adjusting the only un- 

~(2) 

known LEC in the theory, d\ . As discussed in Introduc- 
tion, this is an expected result because the energy depen- 
dence of the experimental total cross section is known to 
be well described by the phase factor and the final-state 
interaction^, which have been taken into account in this 
work, and the overall strength of the cross section can be 
adjusted by the value of d^ . This feature would be the 
same in the NNLO HByPT calculations because an un- 
known constant appears in the contact irNNNN vertex 
and can be adjusted so as to reproduce the experimental 
data though there are many other corrections coming out 
of the pion loop diagrams. 
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